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We investigate the behavior of an A'^-component quantum rotor coupled to a bosonic dissipative 
bath having a sub-Ohmic spectral density J{y>) oc uj" with s < 1. With increasing dissipation 
strength, this system undergoes a quantum phase transition from a delocalized phase to a localized 
phase. We determine the exact critical behavior of this transition in the large-A^ limit. For 1 > 
s > 1/2, we find nontrivial critical behavior corresponding to an interacting renormalization group 
fixed point while we find mean- field behavior for s < 1/2. The results agree with those of the 
corresponding long-range interacting classical model. The quantum-to-classical mapping is therefore 
valid for the sub-Ohmic rotor model. 
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I. INTRODUCTION 

Quantum phase transitions are abrupt changes in the 
ground state properties of a quantum many-particle sys- 
tem that occur when a non-thermal control parameter 
is variedji In analogy to thermal phase transitions, they 
can be classified as either first-order or continuous transi- 
tions. Continuous quantum phase transitions, also called 
quantum-critical points, are characterized by large-scale 
temporal and spatial fluctuations that lead to unconven- 
tional behavior in systems ranging from strongly corre- 
lated electron materials to ultracold quantum gases (for 
reviews see, e.g., Refs. [^-Izl)- 

Impurity quantum phase transitionsS are an interest- 
ing class of quantum phase transitions at which only the 
degrees of freedom of a finite-size (zero-dimensional) sub- 
system become critical at the transition point. The rest 
of the system (the "bath" ) does not undergo a transition. 
Impurity quantum phase transitions can occur, e.g., in 
systems composed of a single quantum spin coupled to 
an infinite fermionic or bosonic bath. Fermionic examples 
include the anisotropic Kondo model- and the pseudogap 
Kondo modelJfi 

The prototypical system involving a bosonic bath is 
the dissipative two-state systemjJAiiS also called the spin- 
boson model, which describes a two-level system coupled 
to a single dissipative bath of harmonic oscillators. Its 
ground-state phase diagram depends on the behavior of 
the bath spectral density J{uj) for small frequencies lo. 
Power-law spectra J(w) oc w'* are of particular interest. 
In the super-Ohmic case (s > 1), the system is in the de- 
localized (disordered) phase for any dissipation strength. 
In contrast, for sub-Ohmic dissipation (0 < s < 1), there 
is a continuous quantum phase transition from a delo- 
calized phase at weak dissipation to a localized (ordered) 
phase at strong dissipation.— In the marginal Ohmic case 
(s = 1), a quantum phase transition exists, too, but it is 
of Kosterlitz-Thouless type i^^i^^ 

The sub-Ohmic spin-boson model has recently at- 
tracted considerable attention in the context of the so- 
called quantum-to-classical mapping. This concept re- 
lates the critical behavior of a quantum phase transition 



in d space dimensions to that of a classical transition in 
d+1 dimensions. The mapping is usually established by 
comparing the order-parameter field theories of the tran- 
sitions: Imaginary time in the quantum problem plays 
the role of the extra dimension in the corresponding clas- 
sical system. In the case of the spin-boson model, the 
classical counterpart is a one-dimensional Ising model 
with long-range interactions that decay as l/r^+'* for 
large distances r. In recent years, the applicability of 
the quantum-to-classical mapping to the sub-Ohmic spin- 
boson model has been controversially discussed after nu- 
merical renormalization group results^^ suggested that 
its critical behavior for s < 1/2 deviates from that of the 
corresponding Ising model. While there is now strong 
evidencei^"— that this conclusion is incorrect and that 
the quantum-to-classical mapping is actually valid, the 
issue appears to be still not fully settled;^^ Moreover, pos- 
sible failures of the quantum-to-classical mapping have 
also been reported for other impurity models with both 
Ising22r— and highe r^^'^^ symmetries; and the precise 
conditions under which it is supposed to hold are not 
resolved. 

In the present paper, we therefore investigate the large- 
N limit of the sub-Ohmic quantum rotor model. Anal- 
ogously to the spin-boson model, this system under- 
goes a quantum phase transition with increasing dissi- 
pation strength from a delocalized phase to a localized 
phase. ^^'^^ We exactly solve the critical properties of this 
transition. Our analysis yields nontrivial critical behav- 
ior corresponding to an interacting renormalization group 
fixed point for 1 > s > 1/2, while we find mean-field be- 
havior for s < 1/2. All critical exponents agree with 
those of the corresponding long-range interacting classi- 
cal model, ■^^ implying that the quantum-to-classical map- 
ping is valid. 

Our paper is organized as follows. We define the sub- 
Ohmic rotor model in Sec. HH In Sec. IIIH we derive its 
partition function; and we solve the self-consistent large- 
N constraint at zero and finite temperatures as well as 
with and without an external field. Section HVl is devoted 
to a discussion of observables and the resulting critical 
behavior. We conclude in Sec. |Vl 
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II. SUB-OHMIC ROTOR MODEL 

A quantum rotor can be understood as a point moving 
on an TV-dimensional hypersphere of radius N^^"^. It can 
be represented by an A^-component vector S satisfying 
= TV. The rotor has a momentum P; the position and 
momentum components fulfiU the usual canonical com- 
mutation relations [SonPp] — iSafs- In the large- limit, 
N ^ oo, the hard constraint = A^ can be replaced 
by one for the thermodynamic average, (S^) = A^, be- 
cause fluctuations of the magnitude of S are suppressed 
by the central limit theorem. The large- A^ quantum ro- 
tor is thus equivalent to the quantum spherical model of 
Ref. [2^ which is given by the Hamiltonian 



1 



2 q2 



hS + f4S^ 



1) 



(1) 



Here, S and P represent the position and momentum of 
one rotor component, ^ is a Lagrange multiplier enforcing 
the constraint (S"^) = 1, and h is an external symmetry- 
breaking field.— 

We now couple (every component of) the rotor to a 
bath of harmonic oscillators.'''^ In the conventional linear- 
coupling form, the Hamiltonian describing the bath and 
its coupling to S reads 



(2) 



with qj, pj, and nij being the position, momentum, and 
mass of the j-th oscillator. The luj are the oscillator 
frequencies and Xj the coupling strengths between the 
oscillators and S. The last term in the bracket is the 
usual counter term which insures that the dissipation is 
invariant under translations in S*.— The coupling between 
the rotor and the bath is completely characterized by the 
spectral density 



2 ^-^ m,w 

j 



(3) 



which we assume to be of power-law form 

J{uj) ^2naujl-'uj'', (0 < cj < . (4) 

Here, a is the dimensionless dissipation strength and uJc 
is a cutoff frequency. We will be interested mostly in the 
case of sub-Ohmic dissipation, < s < 1. 



III. PARTITION FUNCTION AND 
CONSTRAINT EQUATION 

A. Path integral formulation 

We now derive a representation of the partition func- 
tion in terms of an imaginary-time functional integral. 



Because the sub-Ohmic rotor model H = Hs + Hb is 
equivalent to a system of coupled harmonic oscillators 
(with an additional self-consistency condition), this can 
be done following Feynman's path integral approach^! 
with position and momentum eigenstates as basis states. 
After integrating out the momentum variables, we arrive 
at the partition function 



D[SiT)]D[q,iT)] 



The Euclidian action is given by 



As= dT 
Jo 

Ab = 



{S"^ + ojIS^) - hS + fi{S^ - 1) 



(5) 



(6) 
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(7) 



where the dot marks the derivative with respect to imag- 
inary time T, and (3 — 1/T is the inverse temperature. 

The bath action is quadratic in the qj , we can thus ex- 
actly integrate out the bath modes. After a Fourier trans- 
formation from imaginary time r to Matsubara frequency 
ujn, this yields / D[qi{uJn)]exp{-AB) = Z%exp{-AB') 
where Zg is the partition function of the unperturbed 
bath and 



A', 



Z — / / — / 9in 



2m,- w^(a;^ 



■ S{w,,)S{-iu„) . (8) 



The sum over j can be turned into an integral over the 
spectral density J(w). Carrying out this integral gives 



A'b = ir^ac^i-^lu;„1^5(w„)5(-w„ 



(9) 



with the dimensionless coupling constant a = 
27ra cosec(7rs/2). Combining As and A'b yields the ef- 
fective action of the sub-Ohmic rotor model as 



T^/i(a;„)S'(~w„) , 



(10) 



where e = Wq -f 2/i is the renormalized distance from 
quantum criticality. The term in ^5 is subleading 
in the limit ujn 0. It is thus irrelevant for the criti- 
cal behavior at the quantum critical point and has been 
dropped. The theory then needs a cutoff for the Mat- 
subara frequencies which we chose to be Uc- Because 
the effective action is Gaussian, the partition function 
Z = Zg f D[S{ujn)]exp{—Acff) is easily evaluated. We 
find 



Z=ZO exp(/3M)n 



27r 



-,1/2 



Tie- 



K|«), 



X exp 



2 ^ 



h{uJn)h{-UJ„) 



e + aojc 



0J„ 



(11) 
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B. Solving the spherical constraint 



2. T>0 andh^O 



The spherical (large-A^) constraint (S*^) = 1 can be 
easily derived from the free energy F = —T\nZ by 
means of the relation = dF/d^. In the case of a 
time-independent external field h with Fourier compo- 
nents h[ujn) = Sn^oh/T, this yields 



1 



(12) 



We now solve this equation, which gives the renormalized 
distance from criticality, e, as a function of the external 
parameters a, T, and h, in various limiting cases. 



1. T = andh = 

At zero temperature, the sum over the Matsubara fre- 
quencies turns into an integral, and the constraint equa- 
tion reads 



1 



e -I- aujc '^1^'^ 



1 



(13) 



For sub-Ohmic dissipation, s < 1, a solution e > to this 
equation only exists for dissipation strengths a below a 
critical value ac because the integral converges at the 
lower bound even for e = 0. The value of ac defines the 
location of the quantum critical point. Performing the 
integral for e = 0, we find ac = 1/[7t{1 — s)]. As we 
are interested in the critical behavior, we now solve the 
constraint equation for dissipation strengths close to the 
critical one, a < etc. We need to distinguish two cases: 
1 > s > 1/2 and s < 1/2. 

In the first case, the calculation can be performed by 
subtracting the constraint equations at a and at ac from 
each other. After moving the cutoff to oo, the resulting 
integral can be easily evaluated giving 

e = aujcA'^^'-^\ac - ay^^^''^ (s > 1/2) , (14) 

where A = — (1/s) cosec(7r/s). In the case s < 1/2, eq. 
([T3| can be evaluated by a straight Taylor expansion in 
ac — a, resulting in 



acUJcB ^{ac -a) (s < 1/2) , 



(15) 



with B = l/[7r(l - 2s)]. For s < 1/2, the functional de- 
pendence of e on ttc — a thus becomes linear, independent 
of s. As we will see later, this causes the transition to be 
of mean-field type. 

For dissipation strengths above the critical value ac, 
the spherical constraint can only be solved by not trans- 
forming the sum over the Matsubara frequencies in (jl2p 
into the frequency integral in (jl3l) . Instead, the a;„ = 
Fourier component has to be treated separatelyi^ Alter- 
natively, one can explicitly introduce a nonzero average 
for one of the TV order parameter components (see, e.g., 
Ref. 111). Both approaches are equivalent; we will follow 
the first route in the next subsection. 



At small but nonzero temperatures, an approximate 
solution of the spherical constraint (IT2]) can be obtained 
by keeping the a;„ = term in the frequency sum dis- 
crete while representing all other modes in terms of an 
w-integral. This gives 



T 1 



1 



e + auJc ^i^** 



1 



(16) 



We now solve this equation on the disordered side of the 
transition (a < ac), at the critical dissipation strength 
ac, and on the ordered side of the transition (a > ac). 
We again need to distinguish the cases 1 > s > 1/2 and 
s < 1/2. 

In the first case, we subtract the quantum critical (T = 
0, h = 0, a — ac) constraint from ^TE\i . After evaluating 
the emerging integral, the following results are obtained 
in the limit T — > and \a — ad small but fixed. 



a 



■T 



e — A aci^c ^ 
a 

e = eo + 



ac — a 1 — s 



{a > ac, s> 1/2) , (17a) 

(a ^ac, s> 1/2) , (17b) 
T {a<ac, s> 1/2) . (17c) 



Here, eo is the zero-temperature value given in and 
A= —(1/s) cosec(7r/s) as above. For s < 1/2, we expand 
(flBl) in a — ac and find 



a 



a — ac 



■T 



{a > ac, s < 1/2) , (18a) 



e = S-i/2^cC^i/2j,i/2 ^ < 1/2) , (18b) 

e = eo + T (a < ac, s < 1/2) , (18c) 



ac — a 

with eo given in and B = l/[7r(l — 2s)] as above. 

3. T = and h^O 

At zero temperature, but in the presence of an external 
field, the spherical constraint reads 



1 r- , 1 



n Jo e + aujc e 



(19) 



Proceeding in analogy to the last subsection, we deter- 
mine the distance e from criticality in the limit ft, — > 
and \a — ad small but fixed. In the case 1 > s > 1/2, we 
obtain 



1/2 



h 



(a >ac, s> 1/2) ,(20a) 



,l-s^,2s^l/(s+l) 
a s 



e = {A-'acUJc'^h'')^'"^^' (a = a^, s > 1/2) (20b) 



ac — a 1 — s eo 



(a < ac, s > 1/2) ,(20c) 
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where eo is the zero-field value given in ([14]) and A = 
— (1/s) cosec(7r/s) as above. For s < 1/2, the correspond- 
ing results read 



(a > Qfc, s < 1/2) , (21a) 

(a = ac, s < 1/2) , (21b) 
— (a < ac, s < 1/2) , (21c) 



/ \ 1/2 




- OicJ 




1 2 i,2\l/3 


£0 + 







with eo given in (fT5|) and i? = l/[vr(l — 2s)] as above. 



IV. OBSERVABLES AT THE QUANTUM 
PHASE TRANSITION 



After having solved the spherical constraint, we now 
turn to the behavior of observables at the quantum crit- 
ical point. 



A. Magnetization 

The magnetization AI — {S) follows from ((TT|) via M ~ 
-dF/dh = ra(ln Z) jdh. This simply gives 



M = hjt 



(22) 



To find the zero-temperature spontaneous magnetization 
in the ordered phase, we need evaluate (|22p for T = 0, 
a > ac, and /i — >■ 0. Using equations (pOb ) and (PTk). we 
find 



M = \J(pi- ac)/a 



(23) 



for the entire range 1 > s > 0. The order parameter 
exponent (3 thus takes the value 1/2 in the entire s-range. 
For T > 0, e does not vanish even in the limit h ^ {). 
The spontaneous magnetization is therefore identical to 
zero for any nonzero temperature, independent of the 
dissipation strength a. 

The critical magnetization-field curve of the quantum 
phase transition can be determined by analyzing ((22)) for 
T = 0, a = ac, and nonzero h. In the case 1 > s > 1/2, 
inserting ^p) into ^ yields 



M = 



1/(1+'*) 



(s>l/2) (24) 



which implies a critical exponent 5 ~ (1 + s)/(l — s). For 
s < 1/2, we instead get the relation 



2, ,-l!,\l/3 



{s < 1/2) 



(25) 



The critical exponent S thus takes the mean-field value 
of 3. 



B. Susceptibility 

The Matsubara susceptibility can be calculated by tak- 
ing the second derivative of In Z in (ITT]) with respect to 
the Fourier components of the field, yielding 



1 



awl "^Iw,, 



(26) 



We first discuss the static susceptibility Xst = x(0) = 1/e 
in the case 1 > s > 1/2. To find the zero-temperature, 
zero-field susceptibility in the disordered (delocalized) 
phase, a < Oc, we use (jl4p for e, which results in 



Xst 



(s>i/2). (27) 



The susceptibility exponent thus takes the value 7 = 

For dissipation strengths a > ac, the susceptibility di- 
verges in the limit T ^ Q. The temperature dependencies 
follow from substituting (fT7k ) and (|17b ) into Xst = 1/e- 
This yields 



Xst 



.y-i 



a 

Xst = ujl-^a~^A'T~' (a = ac, s > 1/2) . (28b) 

In the ordered (localized) phase, we thus find Curie be- 
havior with an effective moment of — {a — cx^ja in 
agreement with (1^51) . 

The static susceptibility in the case s < 1/2 is obtained 
analogously. Using (fTSj) . the zero-temperature, zero-field 
susceptibility reads 



(a > ac, s > 1/2) , (28a) 



Xst 



a~^uj^^B{pLc - a)-i (s < 1/2) , 



(29) 



implying that the susceptibility exponent takes the mean- 
field value 7 = 1. From (118b ). we obtain the temperature 
dependence of Xst at the critical damping strength. 



Xst = t^c ^"c 



i/2^i/2j.-i/2 (a = ac, s< 1/2) . (30) 



In the ordered phase, the behavior for s < 1/2 is identical 
to that for s > 1/2 given in (pSb). 

We now turn to the dynamic susceptibility. To com- 
pute the retarded susceptibility x('^): "we need to analyt- 
ically continue the Matsubara susceptibility by perform- 
ing a Wick rotation to real frequencies, iaj„ — > w + iO. 
A direct transformation of (pS)) is hampered by the non- 
analytic frequency dependence Iw^l*. We therefore go 
back to a representation of the dynamic term in the sus- 
ceptibility in terms of discrete bath modes [see the action 
(|S])]. As this representation is analytic in w„, the Wick 
rotation can be performed easily. We then carry out the 
integration over the spectral density after the Wick ro- 
tation. The resulting dynamical susceptibility reads 



1 



e + awc '^Y^Y [cos(7rs/2) — i sin(7rs/2)sgn(aj)] 



(31) 
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At quantum criticality (a = ac, T = 0, h = 0), the 
real and imaginary parts of the dynamic susceptibihty 
simphfy to 



cos(7rs/2) 
Rex(w) = TzrT—r- , Inix(a;) 



sin(7rs/2)sgn(aj) 



(32) 



in the entire range 1 > s > 0. Comparing this with the 
temperature dependencies I^Ei^) and (150]) . we note that 
the results for s < 1/2 violate uj/T scaling while those 
for 1 > s > 1/2 are compatible with it. 



C. Correlation time 

To find the inverse correlation time (characteristic en- 
ergy) A = ^j^^, we parameterize the inverse susceptibility 
as e + aa;^^'*|cLi„|'' = e(l + |a;„/A|''). This implies the re- 
lation 



(33) 



The dependence of the inverse correlation time on the 
tuning parameter a at zero temperature and field in the 
case l>s>l/2is obtained by inserting into 
In the disordered phase, a < ac, this gives 

A = tj,A-i/(i-^'(a, -a)i/(i^'^) (s>l/2). (34) 

The correlation-time critical exponent therefore reads 
i^z = 1/(1 — s). Note that this exponent is sometimes 
called just v rather than vz in the literature on impurity 
transitions. We follow the general convention for quan- 
tum phase transitions where v describes the divergence 
of the correlation length while vz that of the correlation 
time. By substituting (|17b) into (1551) . we can also deter- 
mine the dependence of A on temperature at a = ac and 
h = 0. We find A = A^^T. The characteristic energy 
thus scales with T, as expected from naive scaling. 

In the case s < 1/2, the zero-temperature, zero-field 
correlation time in the disordered phase behaves as [using 

m] 



l/s 



(s < 1/2) , (35) 



resulting in the mean-field value vz = l/s for the corre- 
lation time critical exponent. The dependence of A on 
temperature at a = ac and h = follows from ((HId); it 
reads A = 5-i/(2^)a;(2«-i)/(2^Vi/(2s). The characteris- 
tic energy thus scales differently than the temperature, 
in disagreement with naive scaling. 



D. Scaling form of the equation of state 

A scaling form of the equation of state for 1 > s > 1/2 
can be determined by subtracting the quantum critical 
(T = 0, h = 0, a = ac) spherical constraint from the 



general constraint ((12]). After performing the resulting 
integral, we find 



ac 



c2 



T 
e 



(36) 



a 

We substitute e = h/M [from ([H])]; and after some 
lengthy but straight forward algebra, this equation can 
be written in the scaling form 



1/(1-^) 



0, 



(37) 



with X being the scaling function, and r — [a — ac)/a 
being the reduced distance from criticality. This scaling 
form can be used to reproduce the critical exponents /3 = 
1/2, 7 = s/ (1 — s), and 5 = {1 + s) / {1 — s) found above. 

For s < 1/2, the same approach gives a scaling equa- 
tion containing the mean-field exponents /? — 1/2, 7=1, 
and (5 = 3. Moreover, an explicit dependence on the 
cutoff for the Matsubara frequencies remains. 



E. Entropy and specific heat 

Within our path integral approach, thermal properties 
are somewhat harder to calculate than magnetic prop- 
erties because the measure of the path integral explic- 
itly depends on temperature. As the spherical model is 
equivalent to a set of coupled harmonic oscillators, we 
can use the "remarkable formulas" derived by Ford et 
al.,"^- which express the free energy (and internal energy) 
of a quantum oscillator in a heat bath in terms of its sus- 
ceptibility and the free energy (and internal energy) of a 
free oscillator. For our spherical model, they read 



Fs 



Us 



1 



-At- 



du) Ff{u, T) Im 
du) Uf^LJ, T) Im 



d 
d 

duo 



Inx(w) 



Inx(w) 



,(38) 



•(39) 



Here, Ff{uj,T) = T ln[2 sinh(a;/2T)] and Uf{Lo,T) = 
(w/2) coth(cj/2r). The extra —/j, terms stem from the 
spherical constraint. Note that the free energy in ([38)) 
is the difference between the free energy of the coupled 
rotor-bath system and that of the unperturbed bath, 
Fs ^ F - = -T\n{Z/Z%). The same holds true 
for the internal energy. Us ^ U — Ug . 

The frequency derivative of Inx(cj) can be calculated 
from pi|) . giving 



Im 



did 



Inx(w) 



(40) 



esaujt. 



'^w^-i sin(7rs/2) 



[e -I- aujc ''w^ cos(7rs/2)]2 -|- [aujc "uj^ sin(7rs/2)]2 

To calculate the impurity entropy Ss — {Us — Fs)/T, 
we insert (|40l) into ([38]) and ([39| and perform the result- 
ing integral. In the disordered phase, a < ac, the entropy 
behaves as 



Ss = Daujl''T'/eo 



(41) 
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1 > s > 1/2 s < 1/2 

l3 1/2 1/2 

7 s/il-s) 1 

5 (l + s)/(l-s) 3 

uz 1/(1 -s) 1/s 

Tj 2 - s 2 - s 



TABLE I. Critical exponents of the sub-Ohmic quantum rotor 
model. 

in the hmit T — for all s in the sub-Ohmic range 
1 > s > 0. Here, eo is the zero-temperature renor- 
malized distance from criticality given in (|14p . and D 
is an s-dependent constant. Upon approaching criti- 
cality, a ac, the prefactor of the T" power-law di- 
verges, suggesting a weaker temperature dependence at 
criticality. The specific heat can be calculated from 
Cs = T[dSs/dT), it thus behaves as Dsaujl-'T" /eQ. 

We now turn to the critical dissipation strength, a = 
ttc, For 1 > s > 1/2, we find a temperature-independent 
but non-universal (s-dependent) entropy in the limit of 
low temperatures. For s < 1/2, the impurity entropy 
diverges logarithmically as ln(a;o/r) with T — s- 0. In 
the ordered phase, a > ac, we find a logarithmically 
diverging entropy for all s between and 1. 

At first glance, these logarithmic divergencies appear 
to violate the third law of thermodynamics. We empha- 
size, however, that the impurity entropy represents the 
difference between the entropy of the coupled rotor-bath 
system and that of the unperturbed bath. Because the 
bath is infinite, the entropy thus involves an infinite num- 
ber of degrees of freedom and does not have to remain 
finite. Whether the logarithmic divergence occurs only 
in the large- TV limit or also for finite- rotors remains a 
question for the future. 

We note in passing that the entropy of classical spher- 
ical models^ii^ also diverges in the limit T ^ (even 
when measured per degree of freedom). In these models, 
the diverges occurs because the classical description be- 
comes invalid at sufficiently low temperatures. It can be 
cured by going from the classical spherical model to the 
quantum spherical model. This implies that the diverg- 
ing entropy in the ordered phase of the sub-Ohmic rotor 
model is caused by a different mechanism than that in 
the classical spherical model. 

V. CONCLUSIONS 

In summary, we have investigated the quantum crit- 
ical behavior of a large- iV quantum rotor coupled to a 



sub-Ohmic bosonic bath characterized by a power-law 
spectral density J{uj) ~ uj'^ with < s < 1. As this 
model can be solved exactly, it provides a reliable refer- 
ence point for the discussion of more complex and realis- 
tic impurity quantum phase transitions. We find that all 
critical exponents take their mean-field values if the bath 
exponent s is below 1/2. In contrast, for 1 > s > 1/2, 
the exponents display nontrivial, s-dependent values. A 
summary of the exponent values in both cases in shown 
in table HI The exponent rj sticks to its mean-field value 
2 — s in the entire region 1 > s > 0, in agreement with 
renormalization group arguments on the absence of field 
renormalization for long-range interactions.-^ The fact 
that the order parameter exponent f3 is 1/2 in the en- 
tire range 1 > s > is a results of the large- limit; it 
generically takes this value in spherical models. 

Moreover, the behaviors of the dynamic susceptibility 
and inverse correlation time are compatible with uj/T 
scaling for 1 > s > 1/2, while they violate uj/T scaling for 
s < 1/2. We conclude that the quantum phase transition 
of the sub-Ohmic quantum rotor model is controlled by 
an interacting renormalization group fixed point in the 
case l>s>l/2. In contrast, the transition is controlled 
by a noninteracting (Gaussian) fixed point for s < 1/2. 

We now turn to the question of the quantum-to- 
classical mapping. The classical counterpart of the sub- 
Ohmic quantum rotor model is a one-dimensional clas- 
sical Heisenberg chain with long-range interactions that 
decay as l/r^+'' with distance r. The spherical (large- A^) 
version of this model was solved by Joyce;^^ its critical 
exponents are identical to that of the sub-Ohmic quan- 
tum rotor found here. The quantum-to-classical mapping 
is thus valid. 

The properties of our quantum rotor model must be 
contrasted with the behavior of the Bose-Kondo model 
which describes a continuous symmetry quantum spin 
coupled to a bosonic bath. For this system, the quantum- 
to-classical mapping appears to be inapplicable!^ A re- 
lated observation has been made in a Bose-Fermi-Kondo 
model. The main difference between a rotor and a quan- 
tum spin is the presence of the Berry phase term in the 
action of the latter. Our results thus support the conjec- 
ture that this Berry phase term, which is complex and 
has no classical analog, causes the inapplicability of the 
quantum-to-classical mapping. 
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